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INTRODUCTION 


During  the  last  decade,  considerable  research  has  addressed  the  problem 
of  satellite  remote  sensing  of  atmospheric  temperature  profiles.  The 
satellite  observed  radiance  can  be  expressed  through  the  radiative 
transfer  equation  as  a function  of  the  temperature  profile,  the  atmo- 
spheric transmittance,  and  the  amount,  height,  and  optical  properties 
of  clouds  in  the  field  of  view  of  the  instrument.  Measurements  are 
usually  taken  at  a number  of  temperature  sounding  frequencies  where 
the  atmospheric  transmittance  is  assumed  to  be  known. 

Since  the  atmosphere  is  rarely  cloud  free,  it  is  necessary  to 
eliminate  or  account  for  the  effect  of  clouds  on  the  observed  radi- 
ances in  order  to  obtain  the  atmospheric  temperature  profile.  When 
clouds  exist  in  the  field  of  view 

I(v)  - NIcl(J(v)  + (1  - N)  Iclr(v) 

where  I (u)  is  the  measured  radiance,  N is  the  fraction  of  the  field 
of  view  covered  by  clouds,  u is  the  sounding  frequency,  !cld(u)  is 

the  average  radiance  from  the  cloud  covered  portion  of  the  field  of 
view,  and  Ic]r(u)  is  the  average  radiance  from  the  cloud-free  portion 

of  the  field  of  view. 

In  general,  in  the  "solution  of  the  cloud  problem,"  one  attempts  to 
determine  Icjr(u)  from  the  observations.  Two  exceptions  to  this  approach 

are  given  by  Rogers  [1]  and  Chahlne  [2],  Two  basic  approaches  have  been 
suggested  for  dealing  with  clouds.  One  utilizes  measurements  of  a single 
field  of  view;  the  other  employs  multiple  fields  of  view.  Derivations 
and  discussions  of  several  of  these  techniques  can  be  found  in  the  papers 
cited  in  the  bibliography. 

This  report  discusses  a new  multiple  field-of-view  method  which  is 
based  upon  less  restrictive  assumptions  than  Smith's  [3]  model  which  has 
been  used  operationally  by  the  National  Environmental  Satellite  Service 
(NESS),  McMIlllan  et  al . 14]. 

SMITH'S  TWO- FIELD-OF-VIEW  MODEL 

The  following  discussion  is  extracted  from  the  paper  by  Smith  [3]: 

Writing  Eq.  (1)  for  two  different  resolution  elements  subscripted 
by  1 and  2 gives 

M*)  • ",  «cidiN>  ♦ (’  - V Icirl(v>  (2) 
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A solution  for  the  averaqe  clear  column  radiance  can  be  obtained  if 
1^  f N2.  The  fulfillment  of  this  requirement  can  easily  be  determined 

from  the  measured  window  radiances,  I(w),  since  for  this  condition 
I^w)  + I2(w). 

If  the  resolution  elements  are  relatively  small  and  adiacent,  it  is 
reasonable  to  assume  that  = d2^v^  and  *clr^  = *clrl^  = 

Icir2('*)  where  Ic^ r(\>)  Is  the  clear  column  radiance  which  is  assumed  to 

be  constant  throughout  both  resolution  elements.  Under  these  assumptions 
Eqs.  (2)  and  (3)  may  he  solved  to  obtain 


I 


I^v)  - N*I2(v) 

★ 

i - n 


(4) 


★ 

where  N 
Eq.  (4) 


N-|/N2.  N*  may  be  calculated  from  the  window  channel.  From 


* !l(w)  _ Iclr(w) 

N " I2(w)  - I cl r^w^ 


(5) 


The  value  of  Iclr(w)  may  be  computed  from  either  an  in  situ  measurement 
or  from  an  estimate  of  the  surface  temperature. 

The  procedure  for  uslnq  this  model  in  operational  processing  by  the  NESS 
is  described  by  McMillian  et  al . [4]. 


DISCUSSION 

Several  problems  with  the  method  presented  in  the  previous  section  can  be 
identified:  (a)  The  surface  temperature  is  not  well-known,  especially 

over  continental  areas  where  topographic  and  veqetative  differences  can 
exert  considerable  influence;  (b)  Measurement  error  tends  to  be  amplified 

by  the  factor  1/(1  - N );  (c)  The  major  assumption  that  the  average 
radiance  arising  from  the  cloud-covered  (and  also  from  the  cloud-free) 
area  is  the  same  for  both  resolution  elements. 


A less  restrictive  assumption  on  the  behavior  of  the  clear  column 
radiance  will  now  be  considered.  For  notatlonal  curposes,  let  R(u,x,y) 
denote  the  clear  column  radiance  for  the  resolution  element  centered  at 
position  x,  y.  Over  the  geographical  area  of  Interest,  it  will  be 
assumed  that  R(u,x,y)  varies  linearly  in  x and  y,  i.e. , 

R(v,x,y)  = A(v)  + B(v)  x + C(v)  y (6) 

Observe  that  Smith's  assumptions  are  the  degenerate  case  B(\>)  = C(v)  = 0. 
Given  a set  of  R(u,x,y)  the  coefficients  in  Eq.  (6)  can  be  determined 
by  standard  least  squares  curve  fit  techniques.  Unfortunately,  the  set 

R(u,x,y)  is  not  given;  instead  a set  of  measurements,  say  R(v,x,y),  is 
given  from  which  R(v,x,y)  is  to  be  determined. 

Sample  sets  of  measured  radiance  data  for  the  window  channel  of  the 
NOAA-4  VTPR  are  shown  in  Tables  1 and  2.  A rectangular  array  of  seven 
consecutive  scan  lines  and  seven  consecutive  scan  spots  was  chosen  for 
this  analysis. 

The  technique  consists  of  two  steps.  The  first  step  is  to  determine  if 
a given  observation  is  contaminated  by  clouds;  the  second  step  is  to 
synthesize  "clear  column"  radiance  for  those  observations  which  are  cloud 
contaminated. 

An  editing  procedure  is  used  to  determine  those  points  which  are  cloud 
contaminated.  The  following  steps  are  used: 

1.  Points  which  are  obviously  cloud  contaminated  are  determined 
and  flagged. 

2.  Equation  (6)  is  fit  to  the  unflagged  points. 

3.  An  rms  of  deviations  about  the  plane  is  computed.  Only  unflaqged 
points  are  considered. 

4.  Any  point  with  a deviation  below  the  plane  of  more  than  2 times 
the  rms  is  flagged. 

5.  Steps  2 through  4 are  repeated  until  no  new  points  are  flagged. 

Those  points  which  were  not  flagged  are  now  considered  to  be  clear. 
Equation  (6)  Is  now  fit  to  these  remaining  points.  The  value  predicted 
by  the  plane  is  used  to  synthesize  "clear  column"  radiances  for  those 
points  which  were  flagged. 

Various  procedures  are  available  for  step  1.  A procedure  which  works 
well  in  practice  is  based  on  the  observation  by  Chahine  [5]  that 
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(7) 


I 


I 


clr 

clr 


(^i ) - itv,) 

(v2)  - I(v2) 


Xcld(v1 ,Pc)  ' !c1 r(v1  ^ 
Idd(v2’Pc)  " Idr(v2) 


for  a single  cloud  at  altitude  P . (v-|  and  v2  are  two-cloud  sounding 

frequencies  and  I is  the  measured  radiance.  In  practice  v2  = 749  cm-1 

and  v1  = 725  cnT1  have  been  used  for  the  NOAA-4  data.)  For  Pc  < pressure 

of  the  tropopause,  the  left  side  of  Eq.  (7)  is  a monotone  decreasing 
function  of  P with  equality  at  P = Pc«  With  this  in  mind  Pc  = 300  mbs 

is  set  arbitrarily  and  the  observation  at  step  1 is  flagged  if  the  right 
side  of  Eq.  (7)  exceeds  the  left  side.  This  procedure  has  produced  good 
results  in  practice. 

Two  applications  of  this  technique  are  shown  in  Tables  1 and  2.  In  the 
portion  of  the  tables  labeled  "original  data,"  the  underscored  values 
are  those  observations  which  were  flagged  at  step  1.  Underscored  values 
In  the  lower  part  of  the  table  indicate  observations  which  were  flagged 
at  either  step  1 or  4,  and  the  values  which  appear  are  those  which  have 
been  synthesized  to  replace  the  original  observations. 


CONCLUSION 

An  algorithm  for  estimating  clear  column  radiances  has  been  presented 
and  discussed.  The  procedure  is  Iterative  and  uses  the  principles  of 
editing  about  a least  squares  planar  estimate. 

Examples  are  shown  for  data  taken  over  White  Sands  Missile  Ranqe  (WSMR), 
NM,  which  indicate  that  the  procedure  does  a reasonably  goed  job  of 
analyzing  Individual  observations. 
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TABLE  1 

■iOAA-4  VTPR  OBSERVATIONS  FOR  20  FEB  75  CENTERED  OVER  WSMR,  NM 
a.  Original  data 


Spot 


Line 

13 

14 

15 

16 

17 

18 

19 

1 

104.50 

102.75 

99.80 

99.20 

95.40 

95.10 

91.30 

2 

105.95 

105.05 

103.30 

96.00 

96.30 

96.00 

88.70 

3 

106.55 

106.25 

99.20 

101.55 

101.55 

94.85 

86.35 

4 

111.80 

105.65 

103.60 

102.75 

101.25 

92.20 

94.25 

5 

106.80 

107.10 

104.75 

104.20 

99.80 

102.15 

98.90 

6 

108.85 

107.40 

104.45 

107.10 

105.50 

104.75 

98.50 

7 

110.35 

107.10 

105.35 

105.35 

102.15 

98.90 

97.75 

b. 

Data  after 

cloud  editing 

1 

104.50 

102.75 

99.80 

99.20 

95.40 

95.10 

91.30 

2 

105.95 

105.05 

103.30 

100.42 

96.30 

96.00 

94.60 

3 

106.55 

106.25 

103.35 

101.55 

94.85 

95770 

4 

111.80 

105.65 

1(53755 

102.75 

101.25 

98.79 

94715 

5 

106.80 

107.10 

104.75 

104.20 

99.80 

1027T5 

98.90 

6 

108.85 

107.40 

104.45 

107.10 

105.50 

104.75 

98.50 

7 

110.35 

107.10 

105.35 

105.35 

102.15 

98.90 

97.75 

See  text  for  explanation  of  underlined  entries. 


I TABLE  2 

NOAA-4  VTPR  OBSERVATIONS  OVER  WSMR  FOR  9 SEP  75 

a.  Original  data 


Spot 


Line 

7 

8 

9 

— — r.  ~ . . — 
10 

11 

12 

13 

1 

109.65 

121.40 

119.65 

115.25 

100.20 

109.05 

102.55 

2 

106.10 

115.25 

119.65 

105.50 

109.95 

104.65 

102.85 

3 

117.90 

119.65 

119.05 

93.75 

99.90 

100.50 

92.25 

4 

127.00 

125.25 

123.45 

085 

SOS’ 

102.85 

98.45 

5 

120.25 

124.35 

125.25 

lTO? 

90? 

TOO? 

10O0 

6 

119.65 

124.35 

124.95 

1 21 . 70 

07? 

116.40 

117.00 

7 

123.75 

121.70 

123.75 

122.60 

lo377? 

116.40 

118.15 

b. 

Data  after 

cloud  editing 

1 

121.62 

121.40 

119.65 

115.25 

110.04 

109.05 

102.55 

2 

TOO? 

115.25 

119.65 

114.46 

TOO? 

104.65 

102.85 

3 

TT77TO 

119.65 

119.05 

116.09 

113.26 

110.46 

107.64 

4 

127.00 

125.25 

123.45 

TT77TO 

114.89 

112.08 

TOOT 

5 

128.12 

124.35 

125.25 

119.35 

TT67TO 

TOO? 

106.10 

6 

TOOT 

124.35 

124.95 

121.70 

TTOT 

116.40 

117.00 

7 

123.75 

121.70 

123.75 

122.60 

119.73 

116.40 

118.15 

See  text  for  explanation  of  underlined  entries. 
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